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ABSTRACT 


We  consider  the  problem  of  determining  the  shape  of  the 
cross  section  of  a  simply  connected  perfectly  conducting  infinite 
cylinder  from  a  knowledge  of  the  far  field  pattern  for  all  ^mgles 
of  observation  and  small  values  of  the  wave  number.  The  method 
we  propose  relies  heavily  on  conformal  mapping  techniques.  In 
particular  we  show  that  module  the  tremsfinite  dizuneter  each 
Fourier  coefficient  of  the  far  field  pattern  of  the  electric 
field  determines  a  Laurent  coefficient  of  the  conformal  mapping 
taking  the  exterior  of  the  \mit  disk  onto  the  exterior  of  the 
unknown  cross  section.  The  tr2msfinite  diameter  is  determined 
by  changing  the  polarization  of  the  Incoming  wave  and  measuring 
the  far  field  pattern  of  the  resulting  magnetic  field.  Of  par¬ 
ticular  interest  is  the  case  when  only  a  finite  number  of  the 
Fourier  coefficients  of  the  far  field  pattern  are  known,  and  in 
this  situation  we  obtain  error  estimates  by  using  results  on 
coefficient  estimates  for  univalent  functions.  j 
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I.  Introduction. 

In  this  paper  we  shall  survey  and  extend  some  recent 
results  we  have  obtained  on  the  two  dimensional  low  frequency 
inverse  scattering  problem  for  electromagnetic  waves  which  we 
feel  offers  a  considerable  improvement  over  previous  work  in 
the  area.  In  particular  we  are  able  to  (1)  circvimvent  the  need 
for  numerical  analytic  continuation,  (2)  show  how  elementary 
polarization  effects  can  be  utilized,  (3)  reduce  the  original 
infinite  dimensional  problem  to  one  of  finite  dimension  with 
explicit  estimates  given  on  the  error  resulting  from  such  a 
reduction,  and  finally  (4)  show  how  a  priori  geometric  knowledge 
on  the  shape  of  the  scattering  obstacle  can  be  used  to  in^rove 
these  error  estimates.  These  results  are  obtained  by  a  heavy 
reliance  on  the  use  of  conformal  mapping  methods  and  the  theory 
of  univalent  functions  and  hence  are  basically  restricted  to  the 
two  dimensional  inverse  scattering  problem.  However  it  is  hoped 
that  the  insight  thus  gained  in  the  two  dimensional  case  can 
show  the  way  to  further  progress  in  the  case  of  three  dimensions 
(For  the  case  of  acoustic  waves  in  three  dimensions  some  steps 
in  this  direction  can  be  found  in  [2]  and  [4]}.  Our  basic  result 
can  be  described  as  follows.  Consider  the  problem  of  the 
scattering  of  a  plane  time  harmonic  electromagnetic  wave  by  eui 
infinite  perfectly  conducting  cylinder.  Then  from  a  knowledge 
of  the  low  frequency  limit  of  the  first  N  Fourier  coefficients 
of  the  electric  far  field  pattern  polarized  parallel  to  the 
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generator  of  the  cylinder,  together  with  a  single  Fourier  co¬ 
efficient  of  the  magnetic  far  field  pattern  of  a  field  polarized 
perpendicular  to  the  one  originally  considered,  it  is  possible 
to  determine  the  first  N-t-l  Laurent  coefficients  of  the  con¬ 
formal  mapping  which  maps  the  exterior  of  the  unit  disk  onto  the 
exterior  of  the  cross  section  of  the  unknown  cylinder.  Hence  an 
approxiiaatlon  to  the  shape  of  the  cylinder  can  be  obtained,  and 
since  the  mapping  function  is  univalent  error  estimates  can  be 
obtained  by  using  the  Area  Theorem  together  with  special  results 
in  the  geometric  theory  of  functions  (c.f.  [8]).  We  do  not 
discuss  the  error  due  to  the  (finite  dimensional)  problem  of 
computing  the  N+1  Fourier  coefficients.  Such  errors  are 
basically  due  to  the  problem  of  accurately  measuring  the  phase 
and  amplitude  of  the  far  field  pattern  along  with  the  problem  of 
determining  the  Fourier  coefficients  of  the  far  field  pattern 
when  only  part  of  the  far  field  is  known.  For  a  discussion  of 
certain  aspects  of  these  last  two  problems  we  refer  the  reader 
to  [7]  and  C5]  respectively. 

II.  Reduction  of  the  Inverse  Scattering  Problem  to  a  Nonlinear 
Moment  Problem. 

Let  D  be  the  (bounded)  cross  section  of  a  simply  connected 
perfectly  conducting  infinite  cylinder  with  generator  parallel  to 
the  z  axis  such  that  D  contains  the  origin  and  the  boundary  3D 
of  D  is  smooth.  Let  be  the  z  component  of  the  total 


electric  field  due  to  the  scattering  by  the  cylinder  of  an 
incident  time  harmonic  plane  wave  with  electric  vector  given  by 

8  where  the  time  dependence  factor  has  been 

i  Z 

factored  out,  (X/y,z)  are  Cartesian  coordinates,  k  is  the  wave 

number,  and  is  the  unit  vector  in  the  z  direction.  Then 

z 

will  be  a  solution  of  the  scalar  Helmholtz  equation 

A,  E,  +  k^  E,  »  0  (1) 

A  z  z 

2  — 

in  R  \D,  vanishes  on  3D,  and  has  the  asymptotic  behaviour 

E^(t,0)  -  .  i  ei<kr+V4)  ^(6,11)  .  0(r-3/2)  (2) 


where  (r,9)  are  polar  coordinates.  E(6;k)  is  known  as  the  far 
field  pattern  and  our  aim  is  to  determine  the  shape  of  0  from  a 
knowledge  of  E(6;k)  for  small  values  of  the  wave  number  k  and 
0  s  e  ^  2'ir.  In  particular  it  follows  from  Green's  formula  and 
the  asymptotic  behaviour  of  Hankel's  function  (c.f.  C3])  that 


E(9!k) 


-1 


3E^  (p,*) 


expC-ik  p  co8(9-^)]  ds 


(3) 


3D 


where  v  denotes  the  unit  outward  normal  to  3D  and  (Or^)  ere 
the  polar  coordinates  of  a  point  on  3D.  If  we  now  expand  E(6;k) 
in  a  Fourier  series  we  have 


E 


(e;k)  ■  ^  a^ 


(k)  e 


inO 


^  i' 


(4) 
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where 


e„(lc) 


■^u 


3E^ (P»4) 

- 5^5 - 


expC-inO-lkp  cos  (6-^)]  dsd6 


D 

— n  f  3E  (p»^)  -.HnA 

-i  "  -I,; -  .  !"♦  d. 


(5) 


and  J  denotes  Bessel's  function  of  order  n. 
n 


In  order  to  proceed  further  we  need  to  determine  the  low 

3E 

frequency  behaviour  of  z/3v  on  30.  Using  the  method  of 
integral  equations  this  was  rigorously  determined  in  [3]  (see 
also  [€])  with  the  result  that  for  (r,6)  €  30 


3E,  1  3E^  , 

"  Io5Jc  TT  *  0  ^(iogicp^ 


(6) 


where  E^(r,9)  is  the  solution  of  the  following  exterior 
Oirichlet  problem  for  Laplace's  equation: 


E^(r,9)  -  log  i  +  Ej(r,9)  in  R^\D 


A-  E^  -  o  in  R^\D 

2  O 


E^  »  o  on  30 
o 


(7) 


E*  is  bounded  as  r  tends  to  infinity, 
o 
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It  is  possible  to  express  E^(r,d)  in  terms  of  the  (unique) 

conformal  mapping  w^f  (z )  which  maps  the  exterior  of  the  domain 
D  onto  the  exterior  of  the  unit  disk  fi  such  that  at  infinity 
f(z)  has  the  expansion 

f(2)  =  a  2  +  b  +  I  +  ;  a  >  o.  (8) 


Indeed  it  can  be  immediately  verified  that  if  f(z)  is  as  defined 
above ,  then 

EQ(r,0)  -  -log|f(z)|  ;  z  =  re^®.  (9) 

Hence  from  (5),  (6),  (9)  and  the  Taylor  series  expansion  of 
Bessel's  function  we  have  that  for  n*0,  1,  2,... 


an(k) 


i""k“ 

2”n:iogk 


I  ^  log|f  (pe^*)  ds  + 

3D 


0( 


k" 

(logk) 
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Since  (kp)  =  (-1)”  J_jj(kp)  the  seune  expression  (up  to  con¬ 
jugation  and  a  factor  of  plus  or  minus  one)  holds  for  n  s  -1  and 
thus  no  new  information  can  be  gathered  from  the  negatively  in¬ 
dexed  Fovirier  coefficients.  It  can  furthermore  be  easily  verified 
that  the  n“o  term  gives  no  Information  on  the  shape  of  D.  Thus 
if  we  define  (for  n  z  1) 


-  i‘"  2"  n.’  llm  ^®5k 

k-K,  j^n 


(11) 


T7.  wjj.' 


we  arrive  at  the  relation 


Pn  =  I  ^  I09  I  f(pe^'^)  1  p"  e^""^  ds.  (12) 

30 

We  note  that  in  practice  only  a  small  nvimber  of  the  y  can  be 

n 

computed  with  any  degree  of  accuracy  since  for  large  n  small 

errors  in  the  measurement  of  the  a  (k)  will  cause  large  errors 

n 

in  the  due  to  the  factor  of  2*^n.’  in  (11)  (However  see 

Section  III).  Hence  we  shall  assume  that  the  u  are  known  only 

n 

for  n=l,2, . . . ,N+1  where  N  is  a  (small)  positive  integer; 
however  as  discussed  in  the  Introduction  we  shall  asstime  for  the 
purposes  of  exposition  that  these  N+1  numbers  are  known  exactly. 

As  it  stands,  (12)  is  not  too  much  help,  since  what  is 

needed  is  a  relationship  between  and  f”^  (w)  since  then  by 

determining  f  ^ (w)  and  evaluating  it  on  the  unit  circle  |wl“l 
we  can  determine  points  on  3D.  However  from  the  Cauchy-Riemann 
equations  we  have 

^  log  I  f(pe^‘*')  I  “-jI-  arg  f(pe^^)  (13) 

and  hence  from  (12)  we  have 


tSSa-- 


‘n  -  -  M  p-ei"*  ds 

3D 

•-[  ^C£-^w):"d„ 

|w|=l 


-if  ;  cf^cw)]” 

|w|.l 


for  n«l,2 , . . .  ,N+1.  (14)  is  the  nonlinear  "nioment"  problem  of 

the  title  of  this  section. 


The  relationship  (14)  is  unfortxinately  not  sufficient  to 

determine  an  approximation  to  f~^ (w)  for  lw|»l.  Indeed,  as 
will  be  seen  on  the  sequel,  (14)  only  allows  us  to  compute  the 

first  N+1  Laurent  coefficients  of  f (w)  module  the  co¬ 
efficient  a  appearing  in  (8).  The  number  a  has  geometric 

significance  in  that  it  can  be  shown  (c.f.  [8])  that  a”^  is 
the  tramsfinite  diauneter  of  the  domain  D.  In  order  to  determine 
this  quantity  we  need  to  augment  the  relation  (14)  and  this  is 
accomplished  by  trauismitting  a  second  plauie  wave  that  is  polarized 
perpendicular  to  the  one  originally  sent.  In  particular  let 

be  the  z  component  of  the  total  magnetic  field  due  to  the 
scattering  by  the  cylinder  of  an  incident  time  harmonic  plane 

4  If  M 

wave  with  magnetic  vector  given  by  H.  ■  e  6  ,  where  we  are 

1  z 
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using  our  previously  defined  notation  and  the  time  dependent 

factor  has  been  suppressed.  Then  H  will  be  a  solution  of  the 

z 

2  — 

scalar  Helmholtz  equation  in  R  \D  with  vanishing  normal 
derivative  on  3D  emd  having  the  asymptotic  behaviour 


(r.e)  .  .  i  ei(>'«’/4)  ^(9,10  .  (15) 


As  is  the  case  of  the  electric  far  field  pattern  we  assume  that 
the  far  field  pattern  H(0;k)  is  known  for  small  values  of  the 
wave  number  k  and  0  s  6  s  2tt.  We  can  then  write  (Cl]) 


H(0;k) 


j  ^  exp  C-ikpcos  (0-(|))  ] 

3D 


and  if  we  expand  H(0;k)  in  a  Fourier  series 


,•)=)  =  ^  b. 


(k)  e^ 


we  find  that  (Cl]) 


b^(k) 


-i-n  I 


^  e“^"'*’]  ds. 


We  now  note  that  (Cl])  for  (r,9)  c  3D 


H^(r,0)  -  1+ik  H^Cr,©)  +  0  (k^  logk) 
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where  H^(r,0)  is  the  solution  of 


=  rcos0  +  H®(r,0)  in  R^\D 


A-  H  =  0  in  R^\D 
z  o 


av 


=  0  on  30 


is  regular  at  infinity, 

i.e.  in  terms  of  the  conformal  mapping  (8)  we  can  write 


f;o(r,e)  .  i  Re  [fU)  +  ]. 


We  can  now  proceed  exactly  as  we  did  in  the  case  when  the 
vector  was  parallel  to  the  generator  of  the  cylinder.  In 
particular  if  we  define  Yq  and  Yj^  by 


= 


1  ■  b  (k) 
-lim  o 


v,  - 


k.o  -j;r 


we  find  from  (18),  (19),  (21)  that 


Y^  =  Area  of  D 
o 


^1  “  ^  Cl  -  -4  3  f”^Cw)  dw. 


w  »1 


(20) 


(21) 

electric 

(22) 

(23) 


(24) 
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We  now  want  to  use  the  relations  (14),  (23),  (24)  to  determine 

an  approximation  to  (w)  on  lw|*l  (and  hence  an  approximation 
to  3D)  together  with  error  estimates  under  the  assijmption  that  the 
numbers  * '’^N+l,^o'”'^l  known.  This  will  be  the  topic  of 

the  following  section. 

III.  Determination  of  an  Approximation  to  3D  and  Error  Estimates. 

We  now  consider  the  nonlinear  moment  problem  (14)  and  the 
relations  (23),  (24),  and  note  that  there  are  two  main  sources  of 

error  in  trying  to  use  these  relations  to  determine  f”^ (w)  on 
(w|=l.  The  first  is  that  the  are  determined  by  means  of  the 

limiting  relation  (11)  and  since  the  only  known 

approximately  this  can  lead  to  severe  computational  difficulties. 

The  second  source  of  error  is  more  basic  in  that  in  practice  it  is 
only  possible  to  compute  a  finite  number  of  the  ii^,n=l  ,2  , . . .  ,N+1 , 

where  in  general  N  is  small,  and  from  these  N+1  numbers 

together  with  y  and  y.  we  want  to  determine  the  function 
ox  V 

f~^ (w) ,  i.e.  we  clearly  have  insufficient  information  to  determine 

f  ^ (w)  exactly.  Although  we  shall  primarily  be  concerned  with 
this  second  source  of  error,  we  pause  briefly  to  offer  a  few 
comments  on  how  to  reduce  the  error  in  computing  the  To 

this  end  suppose  it  is  known  "a  priori”  that  D  is  contained  in 

a  ball  of  radius  r_  and  that  E„(r,6)  can  be  measured  for 

o  z 
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r  »  r  ,  0  ^  9  s  27r,  Then  expanding  E  (r  . in  a 
o  z  o 

Fourier  series  we  have  that  for  r  i  r^,  0  s  6  s  2it, 


E  (r,9)  .  ^  y 


(25) 


where  (kr)  denotes  a  Hankel  function  of  the  first  kind.  If  we 

n 

now  make  use  of  the  asymptotic  behaviour  of  Hankel 's  fvinction  we 


find  that  the  Fourier  coefficients  a  (k)  of  E  (r  ,0)-e 

n  z  o 


ikr  cose 
o 


are  related  to  the  Fourier  coefficients  E(6;lc)  by 


4  n 


(k) 


(26) 


Therefore  from  (11)  and  (26)  we  have  that  for  nil 


U  =  lim  i  2”"'^  ni  logk 

k"  (kr^) 


=  4irn  r  "  lim  a  (k)  logk 
°  k*»-0  " 


(27) 


(2 ) 

where  (kr)  denotes  Hankel ’s  function  of  the  second  kind^  i.e. 


2n  r  lim  logk 
°  k^O 


n 


2ir 

I 


-ikr  cos 0  . 

e)  -e  ®  3  e^"®  d0  (28) 


for  nil.  In  particular  since 


2ir 

1 

0 


exp  C-ikrcos9+in0  3  d6  *  2iri*'  (kr) 


(29a) 


*  I3?k  Eo(r,e)  +  0( - 3^,), 


(logk)' 


(29b) 


where  E  (r,0)  is  the  solution  of  (7) ,  we  have  that 
o 


2ir 

>”  1 


0 


^770)  d0  ;  n  ^  1 


(30) 


and  (29b) ,  (30)  provide  a  more  practical  method  them  (11)  for 
computing  the  numbers 


We  now  turn  to  the  problem  of  determining  an  approximation 

to  f”'^(w)  on  (w(=l  from  the  relations  (14),  (23),  (24)  under 
the  assumption  that  Yj^  and  the  are  known  exactly  for 

n=l,2, . . . ,N+1.  Prom  (8)  we  have  that  f"^ (w)  has  a  Laurent 
expansion  of  the  form 


f'^ (w) 


(31) 


and  since  3D  is  assumed  to  be  smooth  the  series  (31)  is 
uniformly  convergent  on  |w|»l.  Hence  from  (14),  (31)  we  have 
(computing  the  residue)  that 
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Wl  -  .2^c^ 

2 

P2  -  -2it(c^  +  ) 


(32) 


and  in  general 

u  =  -2irn  a  -i  +  lower  order  coefficients.  (33) 

n  n“i 

Hence,  module  the  transfinite  diameter  a~^,  we  can  determine  the 
Laurent  coefficients  c^  recursively  in  terms  of  the  numbers 

n»l,2,...  .  To  determine  the  transfinite  dieunetei*  we  use 

(24)  and  compute  the  residue  to  obtain 


Yi  -  - 


TT 

a 


(34) 


and  thus  from  (32),  (34)  we  can  determine  a.  Hence  if  we  )cnow 
'  •  •  • ' i^N+1  ^1  determine  . .  .c^,  i.e. 

the  f\inction 


(w) 


(35) 


An  approximation  to  30  can  be  obtained  if  we  evaluate 

fjj^(w)  on  |w|«l  and  the  mean  square  error  in  this  approximation 


'  I- - 


is  given  by 


14 


de 


(36) 


We  now  want  to  obtain  an  estimate  for  the  magnitude  of  this  error. 
From  the  Area  Theorem  in  univalent  function  theory  (C8])  we  have, 
using  (23) ,  that 


y 


o 


(37) 


Hence 


"  iCT  Z. 

n=if+l 


s 


1 


(38) 


An  improvement  on  this  error  estimate  can.  be  obtained  if  it  is 
known  "a  priori"  that  D  is  convex  since  in  this  situation  we 
have  ([8],  p.SO)  that 


2 

I  °n I  ^  an (n+1 ) 


(39) 
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Hence  in  this  case 


E(fi  - 


(n+1)' 


4  ,N+2.2  ^  1 

^  ^  /  7:77:4 

»  5=!i+l 


m 

4  /N+2.2  f 

<N+r>  J 


.  ilJ}±lL__ 

3a^ (N+1)^ 

which  is  a  considerable  inprovement  on  (38)  except  in  the  case 

-2  -1 

when  0  is  a  small  perturbation  of  a  disk  and  we  have  a  = 
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